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ABSTRACT 



The pressure amplitude and phase distribution along 
the interface between a tapered fluid layer and an under- 
lying fast fluid bottom were investigated both theoretically 
and experimentally. 

Two different theoretical models were compared experi- 
mentally: a simple model based on a combination of normal 

modes and ray theory and an exact solution based on the 
method of images. The experiment was conducted at 100 KHz 
with a wedge of silicon oil separated from a large tank of 
fresh water by a thin mylar diaphragm. The simple model 
failed to predict adequately the pressure amplitude and 
phase along the interface. The method of images gave 
accurate predictions. 
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I. 



INTRODUCTION 



As sound from the open ocean propagates up onto the 
continental shelf, lower frequencies become lost to the 
water column' energy is transmitted into the bottom. If 
the bottom has a speed of sound in excess of that of the 
water, then it is predicted [Ref. 1,2] that a highly 
collimated beam will form in the bottom. 

For a channel of constant depth overlying fast bottom, 
normal-mode theory shows that if the frequency of the sound 
is less than the cutoff frequency of the mode energy can 
be transferred into the bottom. This can be seen if the 
normal mode is represented by upward-going and downward- 
going rays in the layer. As the frequency of the signal is 
decreased, the incidence (grazing) angle between the ray 
and the bottom increase. Energy continues to be perfectly 
reflected and remains trapped in the layer with no transfer 
into the bottom, until the incidence angle becomes equal to 
the critical angle. At that frequency the mode stops propa 
gating in the layer and energy is transferred into the 
bottom. This concept can qualitatively be applied to a 
tapered layer over a fast bottom. As a propagating mode 
travels up the wedge, the incidence angle of the equivalent 
ray increases with each bounce until it reaches the critica 
angle, at which point energy will be transmitted into the 
bottom. 
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Theoretical and experimental studies of this phenomenon 
have been carried out at the NPS for three years encom- 
passing three generations of experimental setups. 

Lt. Edwards [Ref. 1] with his thesis adviser. Professor 
A.B. Coppens carried out initial investigations in a large 
tank filled with brine solution over which a thin layer of 
olive oil was carefully added. Since the interface between 
the two fluids was horizontal, a wedge-shaped upper layer 
was created by inserting a board of pressure-release material 
at an angle to the horizontal. The sound field was observed 
in the bottom fluid at various distances and depths. When 
plotted on properly chosen coordinates , the data were 
observed to cluster about a single curve consistent with 
the existence of narrow beam of sound enetering the bottom 
at a relatively shallow inclination. 

In the meantime. Professor Coppens developed a simple 
theoretical model that predicted the shape of the beam 
pattern. 

In an attempt to obtain more precise data, Lt . Netzorg 
[Ref. 2], with his thesis adviser. Professor J.V. Sanders, 
redesigned the apparatus to consist of a wedge of water 
separated from a small tank of brine by a 1-mil Mylar 
diaphragm stretched tightly within a frame. The bottom and 
sides of the tank were lined with absorbing material to 
reduce reflections. The results showed that the entry point 
of the beam was close to where it was predicted by the 
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simple model. Agreement between experiment and theory was 
fair but not good enough to verify the model. 

Subsequent to these first experiments a more refined 
model has been constructed, with silicon oil in the wedge 
and a large water tank for the bottom. But the results 
obtained in this apparatus by Professor Sanders were still 
significantly disturbed by what appeared to be some sort 
of interference. 

Experimental results from the three generations are 
compared to the predictions of the simple model in Table 1-1 
and Fig. 1-1. Significant differences between theory and 
experiment were observed on the lower side of the beam 
pattern. In one extreme example, two peaks of beam were 
observed in a third generation of experiments. Generally, 
the observed angle of depression of the beam is shallower 
and the beam width narrower than predicted by the simple 
model . 

At this point, the research effort was redirected to 
go back to examining the amplitude and phase distribution 
along the bottom; these quantities provide the link between 
the acoustic field in the layer and the formation of the 
beam in the bottom. 

The objective of this research reported herein are to; 

1) Develop a simple, closed form expression for the 
pressure distribution along the bottom that can be 
used to reliably predict the beam pattern to the 
bottom. 
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2) Calculate the pressure (including phase) along the 
bottom predicted by application of the method of 
images and compare with experimental results. 

In this report, two different coordinate systems are used. 

The first is a special coordinate system applicable to the 
simple model; the origin is located at a point on the bottom 
where the lowest mode reaches cutoff. The second is a 
general coordinate system more suitable for the image 
calculations and experiment; the origin is located at the 
apex. These two coordinate systems are shown in Fig. I-2a,b. 
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Table I. Comparison oi' the experiments with the 
simple model 
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BEAM DEPRESSION ANGLE (SIMPLE MODEL) BEAM WIDTH (SIMPLE MODEL) 




Fig. /-2a Coordinate system for the simple model 




Fig. /-2b General coordinate system 
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II. 



THEORY 



A. PRESSURE DISTRIBUTION PREDICTED BY A SIMPLE MODEL 
1 . Basic Concept 

The purpose of this part is to summarize the 
theoretical background necessary to conceptualize the 
propagation of sound in a tapered fluid layer above a fast 
bottom, using a combination of normal modes and ray theory, 
a. Layer of Uniform Depth 

One-dimensional propagation in a fluid layer 
of uniform depth H with speed of sound c^ and density p^, 
overlying a second layer with a faster speed of sound 
c^, and density can be exactly represented by an equation 
of the form 

i (oot-K , x 1 ) 

A 

P = [ A sin(K t z')e n (II-l) 

L n z ' 

n 

n 

where z' is the depth, x' is the horizontal distance, 

K , and K , are the components of the propagation vector 
n ' n 

v t h 

k for the n mode in the z' and x' direction respectively. 
A n is the magnitude of the n*"* 1 mode. For propagating modes, 
the propagation vector k can be represented as in Fig. II-l. 
The angle between k and the horizontal is 6. Snell's law 
gives the relationship between the speeds of sound in the 
layer and the critical angle 9 . 



15 



Fig. II 



As long as 
ever , when 
previously 

k , must be 
z n 

Ref . 3 , 



where 



At cutoff 9 
implies k , 




-1. The propagation vector k and its components 

C 1 

0 = arc cos — (II-2) 

c c 2 

9 < 0 c the energy is totally reflected. How- 
0 > 9 C energy is transmitted through the 
totally-reflecting surface. 

To satisfy the requirements of propagation, 
a constant, real number for each mode and from 

k , 

p_ z 

tan k , H = -(-£)(-— — ) (H-3) 

n I n 



n 



* 



(£-) 2 

c 2 



( II-4 ) 



= 9 c and then the geometry of Fig. II-l 

= — so that 
c 2 
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H 



( II-5) 



k 



Z 



n 




and for a mode with frequency just above its cutoff 
frequency f there is a pressure antinode at the bottom. 
Now by writing Eq. II-l in the form 



P 



, ik , 2 ' -ik , z' 

A z z . , , , . 

I - 1 (e n - e n )e l(a3t-kx ) 

n 



i (art-k , x* +k f z f ) 

. X 2 * 

■ - il l V* 



n 



i (cot-k , x 1 -k , 2 1 ) 

X 2 

- e n n } 



( II-6 ) 



we have two plane waves that propagate with wave number 

— ^ y\ 

k = k ,x ± k iZ . If we assume that the frequency is 
x x 

n n 

slowly lowered, higher modes are cutoff as the frequency 

decreases with k , approaching (n - ^-) 

2 2 ri 

n 

b. A Tapered Layer 

Let us now consider a wedge with angle 6 as 
indicated in Fig. II-2 and II- 3, where X' is the distance 
from 9 = 9 c to the apex, and is the distance of closest 
approach (Turnaround distance) of the indicated ray. In 
this case each bounce off the bottom increases 9 by 28. At 
the point where 9 = 9 c energy begins to penetrate into the 
lower layer. The resulting attenuated ray bounces off the 
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Fig. II-2 . The geometry of a wedge 
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Fig. II-3 . Bounces off the bottom of a wedge 
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top and makes another collision with the bottom, losing 

still more energy. This process continues until the wave 

reaches the "Turnaround distance" x' defined by: 

s 



x' = X' (1 - sin 0 ) 
s o c 



(II-7) 



For the wedge then, the frequency remains constant while 
0 increases with each bounce off the bottom until that depth 
H is reached for which 0=0 where 




n 

From Fig. II-l we obtain at cutoff 



k , = k sin 0 

z ’ c 

n 



( II-9 ) 



where 



k 



2rr 

A 



( 11 - 10 ) 



and substituting Eq. (II-9) into Eq. (II-8) yields 

(n - j) A 
H o 2 sin 0 

c 



( 11 - 11 ) 



Knowing H q allows calculation of the value of x' = where 
energy will begin to be transmitted into the bottom. 
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2 . Pressure Distribution along the Bottom 



A.B. Coppens [Ref. 4] formulated an approximate 
model to predict the characteristics of the transmitted 
beam in the bottom. The expressions for the pressure ampli- 
tude and k^, are: 

A = A q e"* A (11-12) 



k 



x' 





sin 0 - 

( — r~) (H-13) 

1 - — 

1 X' 
o 



where x' , X^, are defined from the geometry of Fig. II-l, 
9^ from Eq. II-2, A q is the pressure amplitude at x' = 0, 
and A defined as 



A 



tan 0 ,2 

£{ b 

3 b2_! 



C r 

[arc 



sin 



2 , 

c +bz 
c (b+z ) 



arc 



sin 



TbT 



] 



1 b 

2 b-1 




[arc 



sin 



2 

c +bz 
c ( 1+z ) 



arc sin c] 



+ 



1 1 
2 b + 1 



c [arc 

i/T? 



sin 



c 

c 



2 



-z 

(1-z) 



arc sin c] } 

(11-14) 



where 



b = 




(11-15) 
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c 



(11-16) 



= cos 8 

c 




z 




(11-17) 



x 1 1 2 

In the case where << tt cos 8 this model may be simpli- 

X o 2 c 

fied further 



P 1 tan6 c 2/2 tX* , 3/2 
P- 6 3 a ' 

A = A e ° (11-18) 



and 



sin 8 - 

( grr) (H-19) 

1 - — 

1 X' 
o 

sin 8 

8 = arc sin ( £) (11-20) 

1 - — 

1 X' 
o 

3 . Pressure Distribution under the Bottom 

In this part we present without detailed derivation 
A . B . Coppens [Ref. 2] work concerning the pressure distri- 
bution under the bottom based on the simple model and 
obtained through the use of the Green's function and the 
method of stationary phase. Two basic assumptions were 
considered for this development. 

a. The plane wave has wavefronts that are parallel 
to the shore line. 
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b. The pressure distribution on the bottom can be 
viewed as an acoustical source which generates 
the beam of sound projected into the bottom. 

The geometry is indicated in Fig. II-4. 

The Green's function formulation for this case can be 
developed through the geometry of Fig. II-4 as follows: 

At the point , z') the instantaneous pressure is given by 

P = P(r',t) = P(r' ) e lwt (11-21) 



where 



? (r ’> = T? I ? s <x "»#^!z'=o dx " (II - 22) 

— CO 

The integrand of Eq . (11-22) is formed from the pressure 

distribution at the bottom 



P s (x") 




A ( x " ) 



-i 



e 




x 
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(II-23) 



and the Green's function for a line source is 



G(r ' t r" ) = i tt H (2) ( |r'-r"| ) 

^ o 1 1 



(11-24) 
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Fig. II-4. 



Plane wave propagation in 



the wedge 
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1 



In Eq. (11-23) the term 



is an a postiori correction 



X' 



' o 

for the cylindrical convergence of sound. In the absence 
of losses into the bottom and if there is no energy loss 
within for the wedge, the pressure amplitude would behave 
as 




(11-25) 



also from the geometry of Fig. II-2 we have 




o 

so that 




(11-26) 



f 



(H-27) 



The term A(x") is the pressure amplitude on the bottom 
and the exponential term is the corresponding phase. 

This pressure distribution is an integral of the 
general form 



oo 

I = / F (x ) e G (x } dx 

— CO 



(11-28) 
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where F(x) varies slowly compared with the exponential 

factor. The value of the integral can be well approached 

with the method of stationary phase. 

Using the above method and assuming the simple 
x 11 1 2 

model for << cos 6 , and for small 3 the following 
o 

expression for the pressure distribution was developed: 



P 




v e 



tan0 c 

s 






V 



(11-29) 



where 



v = ' i 2(1 c ? s e> = (11-30) 

tan 6 tan 0 

c c 

2 . Pressure Distribution along the Bottom Predicted 
by the Method of Images 

Since the simple model did not yield results agreeing 
with the experimental measurements, it was decided to use 
the exact but more cumbersome method of images to attempt 
to predict the pressure distribution on the bottom. 

Let the source be a distance X from the apex of 
the wedge and a depth H below the sea surface as indicated 
in Fig. II-5. 

Define 



A . H 

<j) d = arc sm — 



(II-31) 
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and 



<P Q = 6 - <P d (11-32) 

where 8 is the angle between the sea surface and bottom, 

<j>^ is the angle formed at the apex between the source and 
the sea surface, and <t> Q is the analogous angle for the source 
and the bottom. 

If 4> is the angle formed at the apex between the 
t h 

n image and the bottom then Fig. II-5 reveals 

<+>l = 6 + 4> d = 28 - <p o 

<j > 2 = 3 4> Q + 2 ( B— 4> 0 ) = 28+4) q 

<P 3 = 3<P q + 4(6-<0 o ) = 48~0 O (H-33) 

or more generally 

4> = (n+l)d> + n(S“<f> ) = nS + 4 s for n even 

n o o o 

4) = n<j> + (n+1) (8-4) ) = (n+1) 8 - 4> for n odd 

n o r o o 

(11-34) 
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. 
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This set of equations can be reduced further to one 
expression : 

<J> n = {n + i[l-(-l) n ] }6 + (-l) n 4> 0 = 2 INT [£ji] 3+ (-1) n <J> Q 

( II-35a) 

Also for the symmetric member n' of the paired images 
we obtain 

V = {n ' + HS - (-1) n ' = 2 INT[2^i]B + (-1) n ' cJ> Q 

( II-35b) 

where INT [ ] denotes the largest integer which is equal 

to, or smaller than the argument. 

Considering the geometry of Fig. II-5 and II-6 we 

can derive expressions for the distances r and r . which 

n n 

represents the distance between a receiving point on the 

t h t h 

bottom at distance x from the apex and the n and n' 
image respectively. Also we can derive expressions for 
the angles 0 nQ and 9 , which formed at the receiving point 
between the bottom and the images n and n' respectively. 

These expressions are: 

r n = yX 2 + x 2 - 2Xxcosg n (II-36a) 

r , = l/x 2 + x 2 - 2Xx cosg , ( II-36b) 

n ' V n 
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Fig. II-6. Geometric development of the image' s position 
with respect to the receiving point 
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Let us now define 8 and 9 , as the angles of 

nm n m ^ 

th 

incidence for the m bounce from the bottom for the n 

t h 

and n' image respectively, m = 1 , 2 , 3 , ... (The o bounce 

is the last-one before reaching the receiver.) From the 

geometry of Fig. II-5 and Fig. II-7 we obtain expressions 

for 0 as follows: 
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(11-38) 



The general expression is 



0 = 0 - 2m6 

nm no 



(II-39a) 
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Fig. II-7. Geometric development of reflection angles 
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Also for the symmetric member n' 



of the n fc ^ 



pair we obtain 



0 



n 1 m 



6 , - 2m6 

n ' o 



(II-39b) 



and the maximum number of bounces of the n and n' 
image is 



xM = |(n-l) + i{l - (-l) n+1 } 

j(n’-l) + | 1 -{ (-l) n ’ +1 } 

4*— i 

INT [-ji] = INT [jg~] (11-40) 

If the receiving point is on the bottom, then from Fig. II-8 
we obtain: 

1. The distance from each member of a pair of 

images (n,n') to the receiving point is the 

same , i . e . , r = r , . 

n n' 

2. The number of reflections (from the bottom) 

th 

from the n' image to the receiving point is 

t h 

one more than the n image . 

3. The grazing angle to the bottom is the same 
for each corresponding bottom bounce so that 
the reflection coefficients for corresponding 
bounces are equal, except the last bounce for 
the n image . 
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Fig- II-8. Geometry of symmetric images 
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Consequently R = R , where R and R , are the reflec- 

^ ~nm -n m ~nm ~n m 

th 

tions coefficients for the nr bounce of the n'" and n' 
image respectively. 

t h. 

The reflection coefficients for the n c image are 



given by 
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p l c l 
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nm 
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(11-41) 



where 



~nm 



, , 2.2 2 „ 

1 - (— =-) cos 6 

cl nm 



sin 0 



nm 



for 0 >0 (11-42) 
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(11-43) 
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(11-44) 
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Let us now consider the pressure P(x) at a point x on the 
bottom. The contributions from the paired images are 
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(11-45) 



and 
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th th 

(The n and n' images differ only in that n c image has 

no R „ reflection.) 

-no 

The complex pressure is 
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(11-47) 



A computer program was constructed to calculate this 
complex pressure (amplitude and phase) . This program, 
with the corresponding flow chart, is presented in Appendix A. 

In case of a source so distant that it can be 
assumed to be at r ^ 00 then the geometry of Fig. II-9 
yields the approximation 
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Fig. II-9. The geometry of distance source approximation 
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(11-48) 



6 ~ <f> 

no n 



which simplify the mathematics significantly. Here 



R = R(6 ) = R(4> ) 

-no ~ no - n 



(11-49) 



In this case, the expression for reflection coefficients 
are (see Appendix B) 
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where the right hand sides of Eqs . (11-50) depend 

only on n, then define 
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Also, the geometry of Fig. II-9 gives 



r = r - x cos 4 > (11-52) 

n n 

Substitution of Eq. ((11-50), (11-51), (11-52) into 

Eq. (11-47) yields 



P(x) = 



1 -ikr 

? e 



N 

l e 

n=0 



ikxcos<[) 



n 



INT 



(- 1 ) 



U + 



R(n) 



•}T (n) 



(II-53) 

< 



The last expression for pressure is the approximation for 
a distant source, which is described in detail elsewhere 
[Ref. 4]. 

The advantage of Eq. (11-53) with respect to Eq. (11-47) 
is that because the elimination of product term, the 
calculations are greatly simplified. 
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III. EXPERIMENT 



Initially it was intended that the beam pattern in the 
bottom would be exhaustively investigated, and a larger 
experimental system with silicon oil in the wedge and a 
large water tank for the bottom was designed. But, even 
in the new system, it was observed that the beam pattern 
was being significantly disturbed by unknown interferences 
that caused the apparent angle of the beam to move to 
shallower angles than predicted by the simple theory des- 
cribed in II-A and made the beam seem narrower than predicted. 
This interference was definitely not from the walls or 
bottom of the water tank. It might have been associated 
with the wedge, but whether it was an inherent property of 
wedge propagation omitted from the theory II-A or merely 
an artifact of the experimental setup was the question. 

In an attempt to study this question more thoroughly, 
the experiment was redirected to investigate the amplitude 
and phase of pressure on the bottom of the wedge. 

A. APPARATUS 

The experimental facility consisted of a wedge of sili- 
con oil overlying a large tank of fresh water; the two 
fluids being kept apart by a thin Mylar membrane. A photo- 
graph of the apparatus is shown in Fig. III-l. The sound 
source, positioned at the end of the silicon-oil wedge, 
sent a pulsed signal toward the apex of the wedge. It was 
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possible to position the receiver anywhere within the fresh 
water region, but in this experiment the receiver was mainly 
placed close to the bottom of the wedge. The geometry is 
shown in Fig. III-2. 

1. Water Tank (Fig. III-3a,b,c) 

The water tank, constructed of 3/4 in. thick 
plywood, had the inside surfaces painted with Varathane 
but were otherwise uncoated. It was filled with tap water 
at room temperature. 

2. Wedge (Fig. III-4) 

The silicon oil was contained in a frame whose 
sides were made of 1/2 in. x 2 in. aluminum and whose 
bottom consisted of a thin Mylar sheet. Half-mil thick, 
aluminum-coated Mylar was fastened to the bottom of the 
frame by first taping a sheet of Mylar, metalized side up 
to a flat surface. Then the aluminum frame was placed on 
the Mylar, and the edges outlined by a sign pen. A bead 
of "Royal bond GRIP" contact cement was applied to the 
bottom of the frame and to the marked areas of Mylar. After 
the cement was completely dry (about 15 minutes) , the 
frame was placed on the Mylar and about 50 pounds of weight 
placed on top of the frame. After the contact cement had 
cured for 12 hours, a coat of "Silver print conductive paint" 
was applied to the outside joint between the Mylar and the 
frame. This provides electrical contact between the Mylar 
and the frame. (This was necessary to allow measurement of 
the slope of the wedge as will be explained in a later 
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section.) The excess Mylar was trimmed flush with the 
outside edge of the frame by a razor blade. As a final 
step, a thin coat of the "Silicone rubber sealant" was 
applied to both the outside and inside of the Mylar-frame 
interface to prevent leaks. Because the tension in the 
Mylar tended to bend the walls of the frame inward, a 
3/8 in. diameter aluminum rod was used as a "frame spreader" 
to spread the frame and keep Mylar flat. The length of 
the frame spreader was adjustable. It was placed across 
the shorter dimension of the frame and adjusted until the 
Mylar film was mirror flat. Aluminum arms at each end were 
used to suspend the wedge from the top of the water tank. 

The longer arm was connected to the apex side of the wedge 
and extended in the propagated direction. The slope of the 
wedge was set by raising the end of this arm with wooden 
blocks. The other arm, perpendicular to the propagation 
direction, was connected to the source side of the wedge. 

The horizontal level of the wedge was adjusted by means of 
screws at the ends of this arm. 

3 . Source Transducer 

The source transducer was made by Lt. Netzorg. 
Construction details are discussed in his thesis [Ref. 2]. 
The design of the transducer was quite well done: it has 

sufficient horizontal directivity so that sound is not 
reflected from the side of the frame, there is negligible 
sound propagation in the backward direction, and its active 
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face, 2 cm high by 8 cm wide, is the largest size allowable 
in the wedge. The experimentally measured directivity is 
very close to the directivity predicted for a rectangular 
source (Fig. III-5a,b,c). The theoretical directivity 
is described as follows: 



H ( 6 , <j>) = 



sin ( jkL^sinQ ) 
^■kL^sinO 



sin ( jkl^sintjj) 
ykL 2 sin<J> 



[Ref. 3,5] 



(III-l) 



where 



= Horizontal dimension of an active face 
of a rectangular source. 



~ Vertical dimension of an active face of 
a rectangular source. 

K = — = 4— Wave number 



0 = Horizontal angle measured from a major 

radiation axis. 



<J> = Vertical angle measured from a major 

radiation axis . 



As a matter of convenience, the directivity measurements 
were done in the water instead of silicon oil. The fre- 
quency used in these measurements was adjusted to give the 
same wavelength as produced when the source operates at its 
prescribed frequency in silicon oil. A frequency of 
127.8 kHz in fresh water (c 2 = 1503.8 m/sec) produces the 
same wavelength as 100 kHz in a silicon oil (c^ = 1176.5 m/sec) . 
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For more accurate comparison between the experimental 
results and the method of images, the vertical directivity 
was employed in the computer model "WEDGE 0" (APPENDIX-A) . 

The source transducer was mounted on a support that 
allowed it to be located anywhere in the upper fluid layer. 
Vertical adjustments were made by means of a micrometer with 
1/2 in. of travel. 

4. Receiver and Attachment (Fig. III-6a,b) 

A LC-10 hydrophone, manufactured by Celesco Indus- 
tries, was used as the receiver. It was possible to locate 
the receiver anywhere in the fresh water. Since the main 
purpose of the experiment was to investigate the sound 
pressure distribution on the bottom, the hydrophone was 
usually located as close to the bottom as possible. The 
attachment mounted on the frame of the wedge as shown in 
Fig. III-6a . 

At the beginning of this thesis, some measurements 
were made of the beam pattern in the lower fluid. Then a 
different attachment was used (Fig. III-6b) which was mounted 
on the large arm. 

In both cases, the depth of they hydrophone was 
adjustable. In all cases, measurements were made in the 
vertical plane that contained the acoustic axis of the source 
transducer . 

5. Slope-measurement Assembly (Fig. III-7a,b,c) 

The slope measurement assembly consisted of a pin, 
micrometer, and attachment. Fine vertical adjustments of 
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the pin were made by means of the micrometer with 1/2-in. 
of travel. The attachment was mounted on the wedge frame 
and could be located anywhere in the wedge. The pin was 
electrically insulated from the frame, so that when an 
Ohmmeter was attached between the frame and the pin, the 
precise depth at which the pin contacted the metalized 
Mylar could be determined. 

6 . Electrical Setup 

Figure III-8 is a block diagram of the electrical 
setup which consists of a driving system, a receiving system, 
and a measuring system. The transducer was driven by a 
General Radio Oscillator (Type 1310) and Tone Burst Generator 
(Type 1396-A) , polarized by a 300 Volts D.C. biasing-voltage 
[Ref. 2]. The driving frequency (100 kHz) was measured with 
a Hewlett Packard electronic frequency counter (Type 5233L) 
and the driving waveform displayed on one channel of a 
Tektronic Dual-Beam Oscilloscope (Type 565) . The received 
signal was amplified (40 dB) by a Hewlett Packard Amplifier 
(Type 465-A) , and band passed (40 kHz - 200 kHz) through a 
Spencer-Kennedy Laboratories Band Pass Filter (Serial 1683 
Model 302). The received signal was displayed on the other 
channel of the Tektronix Dual-Beam Oscilloscope (Type 565) 
and the amplitude was measured on this display. The received 
signal also went to a Hewlett Packard Oscilloscope (Model 
120B) where its phase was compared to that of the driving 
signal . 
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7 . Auxiliary Apparatus 



To measure the speeds of sound and densities of 
fluids, other pieces of equipment were used. 

Transducer 

Two Hydrophones LC5-2 
Analytical Balance 
ASTM Vessel 

B. PROCEDURE 

Before beginning the experiments, the speeds of sound 
and densities of the fluids (a silicon oil for the upper 
layer, fresh water for the bottom) and the slope of the 
bottom was measured with high accuracy. 

1. Sound Speed Measurements (Fig. III-9) 

The speeds of sound were measured in a pan which 

« 

had dimensions 10 3/4 in. x 7.0 in. x 1 3/8 in. About 1500 cc 
of a test fluid was carefully poured into a pan to prevent 
the entrapment of air bubbles. Two LC5-2 receivers and a 
transducer were used to measure the speeds of sound. The 
distance between two receivers was set at exactly 20.0 cm. 

And the transducer was driven by a 150 kHz pulsed signal. 

The time of transit for one receiver to the other was measured 
using the time-delay feature of the dual-beam oscilloscope. 

The speed of sound was then calculated by dividing the 
distance by the transit time. 

2 . Density Measurments 

The densities of the two fluids were found by first 
weighing a 100 ml, ASTM vessel on an analytical balance 
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readable to one-tenth mg. The vessel was then filled to 
the 100 ml mark and weighed again. By subtracting the first 
weight from the second weight, the weight of 100 ml of the 
test fluid was determined. The weight of the test fluid in 
grams divided by the 100 ml volume and, the quotient times 
1000, gave the density of the test fluid in kg/m^. 

3 . Measurement of the Slope Angle 

The approximate slope was set by raising the end 
of the long arm a precalculated height with the use of 
wooden blocks. Then 3500 cc of silicon oil were slowly 
poured into the wedge, and the frame-spreader was adjusted 
to keep the bottom flat. The water level in the tank was 
then adjusted to keep the apex line (shore line) straight. 

The source transducer was set at about middle depth on the 
wedge centerline. The slope angle was accurately measured 
with the slope-measurement assembly described in the appara- 
tus section. Because the top and bottom of the frame of 
the wedge are parallel, if the Mylar sheet is flat the 
distance from a reference level to the bottom should be the 
same at all points. The distance from a reference level to 
the surface and to the bottom were measured at several points 
on the wedge centerline, and the data were plotted. The 
slope was calculated from the plot as shown in the example 
illustration in Fig. III-10. 

4 . Measurement of Pressure Amplitude along the Bottom 

The electrical setup of the measurements was shown 

in Fig. III-8. The oscillator output, set at 100 kHz, was 
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fed into the tone burst generator. A pulsed signal of 
36 cycles of the 100 kHz input with a closed-gate duration 
time of 10 msec (repetition rate was about 100 Hz) went 
through the transducer biasing network and then to the 
source transducer. The driving signal was displayed on one 
channel of the dual-beam oscilloscope and adjusted for a 
peak-to-peak voltage of 10 Volts. The receiving signal was 
displayed on the other channel of the same dual-beam oscillo- 
scope. This was very convenient because the receiving sig- 
nal was always compared with the driving signal on the same 
oscilloscope. Reflections of the receiving signal were not 
observed. 

5 . Measurement of Pressure Phase along the Bottom 

The phase of the receiving pressure was compared to 
that of the driving signal on an oscilloscope. Fig. III-8 
also shows the electrical setup for these measurements. The 
signal from the oscillator was directly fed into the x-axis 
of the oscilloscope. The receiving signal was amplified 
and filtered, then fed into the y-axis of the same oscillo- 
scope. The signal on the x-axis was continuous with a con- 
stant amplitude (adjusted to 10 Volts peak-to-peak voltage) . 
On the other hand, the signal applied to the y-axis was 
pulsed and had an amplitude that depended on the receiver 
position. The pattern observed on the oscilloscope was an 
ellipse crossed by a line on the x-axis (like "9" of the 
Greek alphabet) . The line occurred during the time the 
received pressure was zero, i.e., between pulses. When the 
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signal, generated by the oscillator, was sent directly to 
the transducer, bypassing the tone burst generator, the 
pattern on the oscilloscope was the same ellipse but without 
the line on the x-axis. This indicates that the pulse is 
sufficiently long to allow accurate phase determinations. 

(It also showed that reflections in the tank were negligible.) 
But measurements were still difficult, because the amplitude 
of the receiving signal was not constant and the phase 
changed very quickly with changes in the receiver position. 

The phase measurements were done independent of the ampli- 
tude measurements. The position of the receiver was recorded 
when the phase was 0°, 90°, 130° and 270°. (Corresponding 
to a straight line with positive slope, a symmetric ellipse 
and a straight line of negative slope.) 
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Fig. /Jl- 2 Geometry of the system 
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Fig.ii /-3a Water tank and wedge (Top view) 
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Fig.JJJ-3b Water tank and wedge (Side view l) 
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Fig.//J-3c Water* tank (Side view 2) 
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Fig. Ul~5a Directivity of the source transducer 




Fig. ///-5b Horizontal directivity pattern H(9, 0) 
and boundary of the wedge 
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Fig. I! /-5c Vertical directivity pattern H{o,<p) 
and boundary of the wedge 
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SCALE 3/10 



Fig J//-6a Receiver and. attachment (l) 
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Fig.ff/-6b Receiver and attachment (2) 
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Fig.n/-7b Slope-measurement assembly and wedge 

(Side view) 
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Tig. ///-7c Slope-measurement assembly and wedge 

(Side view) 
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Fig, 1 11-8 Electrical Setup 
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DISTANCE FROM THE SOURCE (cm) 
Fig.jjj-io Measurement of the slope (CASE 2) 



IV. RESULTS AND COMPARISON 



Measurements were conducted under two different 
conditions : 

Common Parameters 



Speed of sound in the wedge 
(silicon) 

Speed of sound under the bottom 
(fresh water) 

Density of a silicon oil 

Density of a fresh water 

Frequency 

Pulse length 

Closed gate duration 

Pulse repetition rate 



c 1 = 1176.5 m/sec 



c 2 = 1503.8 m/sec 



= 1015.0 kg/m' 
p, = 995.0 kg/irf 



f = 100 kHz 



= 36 wavelengths 

t =10 msec 
c 

f = 100 Hz 
P 



CASE 1 



Slope angle 

Source distance from the apex 
Source depth 



6 = 0.047 rad 

X = 7 0 . 0 cm 
H = 1 . 6 4 cm 



CASE 2 

Slope angle 

Source distance from the apex 
Source depth 



3 — 0.042 r ad 

X = 75.5 cm 
H = 1.57 cm 



The speeds of sound and the densities of the fluids were 
measured at 18.5°C. The two cases differed mainly in the 
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slope angle: for the first case 3 = 2.7° and for the 

second 3 = 2.4°. These are fairly realistic values: 
according to Sverdrup-Johnson ' s - Fleming [Ref. 6], con- 
tinental slopes off mountainous coasts have, on the 
average, a slope of about 3.5°, whereas off coasts with 
wide, well-drained coastal plains, the slope is about 2°. 
However, there are many different estimates based on differ- 
ent measurements at different locations. For extreme 
examples, estimates of the average continental slopes range 
from about 5° to about 0.1° [Ref. 7], 

The distance from the apex to the source was limited 
to less than about 6 times the lowest possible cutoff 
distance x q (x q = 9.92 cm for CASE 1, x q = 11.34 for CASE 2). 

A. EXPERIMENT VS. SIMPLE MODEL 
1. Pressure 

Figures IV-la,b show the experimentally-obtained 
pressure distribution along the bottom. In both cases 
there was a dominant pressure peak near the lowest cutoff 
point but shifted slightly in the direction away from the 
apex. To compare the simple model with the experiment, 
data and the theoretical curve were plotted on the same 
graph for the interval 0 <_ x <_ 2x q (Fig. IV-2a,b) . There 
are significant differences between the experimental results 
and the simple model. The theoretical curves decay quickly, 
and the model predicts that there should be no energy 
transmitted into the bottom, outside the region x g x <_ X o 
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(the small amount observed in the simple model results 
from the approximations) . On the other hand, the experi- 
mental curves have a long tail with energy transmitted 
continuously into the bottom all the way to the apex, and 
a pressure peak shifted significantly, so that energy is 
leaked into the bottom before the cutoff depth, x > x Q . 

2 • Phase 

The phase distribution was measured in the region 
0 < x < 6x q for CASE 2. The rate of change in the measured 
phase was almost constant but with small irregular fluctua- 
tions. Figure IV- 3 compares the experimental results within 
the region 0 _< x <_ 2x q with the results of the simple model 
calculated from Eq. (11-13) 



kx' 




• 2 o 

sin 0 

c 




Phase = -k ,x' 
x 




sin 




8 

c 

2 



(IV-1) 



For comparison with the experiment, the phase was restricted 
to the interval -it to +ir. 

Figure IV- 4 shows the apparent angle of incidence 

sin 6 sin 0 

c c 

0 = arc sin ( — r) = arc sin ( ) 

2 ^ x x 

^ X o (IV -20) 
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which is calculated from the phase (-k ,x') as discussed 

X 

in Section II-A2. The graph shows that comparison between 
the simple model and the experiment is better accomplished 
by studying 0. The rate of change in the phase is almost 
constant when x >_ x q , but the incident angle 0 changes 
significantly in this region. The incident angle 0 increases 
with decreasing distance x. At x = x q the incident angle 
0 becomes the critical angle and Eq. (IV-1) shows that 
-k ,x' = 0. Subsequently the incident angle increases more 

X 

rapidly as x decreases from x to x . At the turnaround 

os 

point x = x the angle 0 becomes 90° and again -k ,x' = 0. 
Thus, energy can be transmitted into the bottom only for 
x g £ x <_ x q according to the model. Precise experimental 
investigation of the incident angle was difficult because 
the rate of change in the phase was very small. Significant 
differences between experiment and the simple model were 
observed. From the experiment, there is no turnaround 
point. The experimental pressure distribution has a long 
tail which means energy can transmit into the bottom all 
the way to the apex. And the slopes of the phase curves 
(Fig. IV-3) show significant differences between experiment 
and the simple model. 

B. EXPERIMENT VS. METHOD OF IMAGES 

The computer program "WEDGE 0" was designed for a com- 
parison between experiment and the method of images. The 
program has the flexibility of an adjustable source distance 
and inclusion of the vertical directivity of the source. 
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1 . 



Pressure 



Before beginning the comparisons, the computed 
pressure distributions were observed with several combina- 
tions of distance (x = 100 m) or near source (x = 1 m) and 
with directional or omni-directional source. In all cases 
the major pressure peaks and tails were similar; the 
distance-source case showed slightly smoother pressure 
distributions and was little effected by source directivity. 

On the other hand, directivity becomes important if the 
source was near. As a result, computations were done with 
the actual source distance (x = 70.0 cm for CASE 1, 
x = 75.5 cm for CASE 2) and source directivity described 
in Section III-A3, and these results were compared with 
the experiment. 

Figures IV-5a,b show the comparison between experi- 
ment and calculation. Both pressure amplitudes were normalized 
at the maximum pressure peak. In both cases, there are some 
quantitative differences, but good qualitative agreement ; 
most of the local maxima and minima occur at the proper 
positions and with the approximately correct amplitudes. 

Figures IV-6a,b are enlarged graphs of Fig. IV-5a,b 
for the region 0 <_ x <_ 2x q . From these graphs, it is 
obvious that the tails of major peaks of the experiment 
and theory agreed very well. 
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2 . 



Phase 



Theoretical calculations were also done with the 
computer program "WEDGEO" for a distant source. The 
results were compared with experimental data (Fig. IV- 7, 

CASE 2 only) with both phases normalized at x = X . In 
the region X Q <_ x <_ 2X q , these phases were almost identical, 
but some differences were observed for x < X q ; the experi- 
mentally observed phase changed somewhat more rapidly 
than that of the theoretical prediction, the difference 
increasing as the apex was approached. There are some 
possible explanations for this. The depth of the wedge 
at x = X Q was 4.7 mm, and most differences were observed 
where the depth was 3.0 mm or less (0 <_ x <_ 6 cm) . In 
this region, the pressure amplitude was very small 
(Fig. IV- 2a) , so that measurement of the phase was quite 
difficult. Also, very close to the apex the surface of 
the silicon oil might be affected by surface tension and 
not remain ideally flat. 

The incident angles 0 in the upper layer were approxi- 
mated as follows; from Fig. II-l, the angle of incidence 
is, 

k , 

6 = arc cos (■■■: * ■■■ ) 

k 



where 



k 



(jo 2irf 
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■' 





x_ for which k ,X' have the 
2 x ' 



Two adjacent points x^ , 
same values in the interval -tt to +ir yield a change in 
phase of 2tt, 



-k ,X' + k ,X' 

x 2 2 ^ 1 



= 2tt 



k = 2 tt 
x' x| - X 2 



2ir 

Ax 



and 



0 = arc cos (f^) 



(IV- 2) 



Figures IV-3a,b show the results for the experiment 
compared to the results obtained from the method of 
images and from the simple model. There is no turnaround 
point (for which we would have 0 = 90°). In both cases 
the incident angle equals the critical angle a little 
further away from the apex than x/X q = 1, which is consis- 
tent with the measured behavior of the amplitude. For 
x > X Q , most of these incident angles were less than the 
critical angle 0 c , (This does not imply that there is no 
energy transmitted into the bottom: these incident angles 

result from the phase-coherent interferences of a large 
number of images, the fields of which have angles of 
incidence on the bottom encompassing a range of values, 
some of which exceed © c . This shows a significant differ- 
ence between the method of images and the simple model for 
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which energy can not be transmitted into the bottom when 
the incident angle is less than the critical angle.) 

The phase as predicted by the method of images and 
the experimental phase is nearly linear in x. The ampli- 
tude varies more slowly, but perhaps not slowly enough 
for the method of stationary phase to be applicable. 



C. ATTEMPTED RECONCILIATION 



The predictions of the simple model do not agree well 
with the experimental measurements. However, investigations 
of the pressure field in the bottom would be facilitated 
if a simple analytical expression for the pressure distri- 
bution along the wedge-bottom interface could be formulated. 

Consequently, the purpose of this section is to com- 
pare the simple model with the method of images for those 
cases where the source is very distant. The amplitude of 
the pressure distribution predicted by the simple model 
has the form 



A = 



A o e 



Pi k i tan0 c k 2 x' k 4 

-<P7> • k 3<£- ) 

2 o 



( IV- 3 ) 



2/1 



where the predicted values are k^ = - 1, k^ = — 



and 



k^ = . Changing p^, P 9 , c^, c^ and b it was found that 

the pressure amplitude depended only on the ratios — , 
tan 0 p 2 



Next comparing the predictions of the method of images 
with the simple model also showed that the value k^ = j 



12 



gave the proper spatial behavior, to reasonable accuracy. 

This was accomplished by taking In of both sides of 

A 

Eq. (IV-3) and plotting then In vs Q n log-log paper. 

o 

Obviously the slope of the equation 



ln{ln(-^) } 



P 1 k l 

in! (A, 1. 



tan 0 k , 

( — g— E ) -k 3 } + k 4 In (2-,) 



( IV- 4 ) 



gives k^ . 



By the same technique the exponents k, and k„ 
tan 6 



Pi 

of — 
Po 



and 



and the factor k_ were examined but in these 

3 3 

cases the results for k^, , k^ were not constants; but 

P i c i x ' 

functions apparently dependent on — , — , 3, rrr- 

p 2 c 2 A o 

We can therefore conclude that the simple model is 
not generally valid for arbitrary values of p^, , 

X 1 

> 3 and However, it is useful to note that by assuming 

1 ° 

k^ = j the simple model predictions are brought into good 



agreement with the experimental measurements for the values 



3 

= 1.015 gr/cm 

3 

p 2 = 0.995 gr/cm 

c^ = 1173.5 m/sec 

C 2 = 1503.8 m/sec 
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38.706° 



e 

c 

S = 0.0476 rad 



After the discovery of this agreement an effort was made 
to specify the ranges of variables in which the amplitude 
prediction of the simple model is valid with = j. 

The simple model requires 



x 1 1 2 

T- « I cos e c 

o 



( IV- 5 ) 



Assuming that the amplitude distribution is unimportant 

x 1 A 1 

for those values tz-t such that - — << — , this ineauality must 

X A e 2 

o o 

be achieved within the restriction (IV-5) . 



p. tan 9 
1 c 



_A_ 

A 



= e 



8/2 



2 q 
cos 9 



<< 



-1 



(iv-6: 



or 



1 2 M 1 

8 << sin 9 (1 - sin 9 ) — 

8/2 C c P 2 



( IV- 7 ) 



This inequality was tested on the interval 0.002 <_ S £ 0.005 



for 9 =60° and — 

c p. 



0.5 or 1.0 there was uniformly good 



agreement between the simple model and the method of 



images for 



= 0.5, but for 



= 1.0 there was strong 



disagreement between the two models. 
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Therefore the range of validity of the simple model 
with ^3 = ^ appears badly limited. Further, since the 
phase of the simple model does not agree with the experi- 
mental measurements and the predictions of the method of 
images (as already has been mentioned in part IV-A2) , 
we can conclude finally that this model is very weak and 
its application very restricted. 
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NORMALIZED DISTANCE A'/Xo 
Fig. /\/- la Experimental pressure distribution CASE 
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NORMALIZED DISTANCE .\/\' o 
Fig.iV-lb Experimental pressure distribution CASE 2 
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NORMALIZED DISTANCE X/X 0 

Fig.H-2a Comparison of pressure amplitude CASE 
(Experiment VS. simple model) 
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Fig./V-Sb Comparison of pressure amplitude CASE 2 
(Experiment VS. simple model) 




(•pej) asVHtl 



80 



Fig.iV-3 Comparison of phase CASE 2 

(Experiment VS, simple model) 
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A's/Ao 1.0 

NORMALIZED DISTANCE A/A'o 

Fig .IV - 4 Theoretical phase and incident angle (Simple model) 
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NORMALIZED DISTANCE a/Xo 

Fig./V-5a Comparison of pressure amplitude CASE 
(Experiment VS. method of images) 







VO 



CO 



(N 



o 



d aanirHMv sHnssstfd cEuznvwHON 



83 



NOKMALIZED DISTANCE X/Xo 

Fig./\-5b Comparison of pressure amplitude GAGE 2 
(Experiment VS. method of imafes) 
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NORMALIZED DISTANCE X/X 0 

Fig.iV“6a Comparison of pressure amplitude CASE 
(Experiment VS. method of images) 




d aaaindwv annsssid: aaznvwHOM 
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NORMALIZED DISTANCE v/x 0 

Eig.lV-6b Comparison of pressure ainplitude CASE 2 
(Experiment VS. method of images) 




i 
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Fig./V-7 Comparison of phase CASE 2 METHOD OF IMAGES 

(Experiment VS. method of images) 
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NORMALIZED DISTANCE A/Xo 

Fig. /V-8b Comparison of incident angle CASE 2 
(Experiment VS. theories) 



V. CONCLUSIONS 



Two main conclusions should be emphasized. 

(1) After a careful examination of the simple model 
it is found to fail to describe adequately the pressure 
along the bottom and it is therefore not useful for future 
efforts in describing the pressure radiated into the 
bottom. 

(2) On the other hand the method of images gives a 
good prediction of pressure amplitude and phase distribu- 
tions along the interface and therefore offers us hope in 
attempting to predict the properties of the beam of sound 
in the bottom. 
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APPENDIX A 



A computer program for the calculation of the pressure 
amplitude and phase distribution along the bottom by the 
method of images. 
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FLOW CHART OF THE SUBROUTINE "FRES" 
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c 

c 

c 

c 

c 



WEDGE 0 

*** MAIN PROGRAM *** 

r. 
* 
* 
¥ 
¥ 



TO OBTAIN 1 Hl£ PRESSURE 
A WEDGE-SHAPED FLUID 



¥ THIS PROGRAM EMPLOYS THE METHOD OF IMAGES 
*AND PHASE DISTRIBUTION ALONG I HE BOTTOM OF 
*LAYER, OVERLYING A FAST BOTTOM ♦ 

¥ 

¥ ASSUMPTION 

¥ PLANE DAVE REFLECTION COEFIC IENTS . * 



* NOTATION 

I/O OUT-OUTPUT 
IN = 1 N P U r 
CHARACTER 

V VARIABLE 
P -PARAMETER 
C -CONSTANT 

R -REAL NUMBER 
I = INTEGER 
Z -COMPLEX NUMBER 



SYMBOL 



MEANING 



I/O CHAR AC H R 



PA PRESSURE AMPLITUDE ON THE BOTTOM OUT 

PH PRESSURE PHASE ON THE DOTTOM OUT 

Cl SOUND SPEED IN MEDIUM 1 (WEDGE) IN 

C2 SOUND SPEED IN MEDIUM 2 (BOTTOM) IN 

ROl DENSITY OF MEDIUM 1 IN 

R02 DENSITY OF MEDIUM 2 IN • 

BETA SLOPE ANGLE OF BOTTOM IN 

X SOURCE DISTANCE FROM TH F APEX IN METERS IN 

XX DISTANCE FROM THE APEX IN METERS OUT 

AT WHICH PA AND PH ARE CALCULATED 
F DRIVING FREQUENCY IN 

H SOURCE DEPTH IN METERS IN 

ANGLC CRITICAL ANGLE 

HO LOWEST POSSIBLE MODE CUT OFF DEPTH IN METERS 

XO LOWEST POSSIBLE MODE CUT OFF DISTANCE FROM THE 

APEX IN METERS 

XN1 DISTANCE FROM APEX AT WHICH CALCULATION STARTS IN 
NORMALIZED DY DIVIDING BY XO 

XN2 DISTANCE FROM APEX AT WHICH CALCULATION STOPS IN 

NORMALIZED BY DIVIDING DY XO 
DISTX1 DISTANCE XN1 IN MEIERS 

DISTX2 DISTANCE XN2 IN METERS 

DX RANGE OF CALCULATION IN METERS 

N NUMBER OF POINTS FOR WHICH PA AND PH IN 

CALCULATED 

PA I 3.1415 

WN WAVE NUMBER 

C21 SOUND SPEED RATIO -C2/C1 

RC21 ACOUSTIC IMPEDANCE RATH) 

ANGLO ANGLE FORMED BY SOURCE > AFEX AND BOTTOM 
DISTX TEMPORARY FOR XX(I)/X 
RAMP TEMPORARY FOR PA 

PHAS TEMPORARY FOR PH 



92 



nnnn o n n n nonnnn n n o non non nnnnnn 



< (PROGRAM) ) 

DEFINE CHARACTERS OF VARIABLES rPARAME TERS AND CONSTANTS 
IMPLICIT REALMS (A-HrG-Z) f INTEGER ( 1-N) 

DIMENSION F*A< lOi ) rPH< 101 ) rXX( 101 ) 

COMMON C21 r RC21 r BE T A 7 UN r PAI f ANGLO 7 X f H »F 

READ IN DATA 

READ ( 5 f 1 00 ) Cl fC2fR01 »R02»BETA 
READ (5, 500) XN1 >XN2, X 7 H 7 F 
READ (5 r 600) N 

CALCULATE AND WRITE CONSTANTS 

WRITE ( 6 ? 700 ) C1fC2fR01fR02»XN1fXN2fXfHfBETAfF 
ANGLC — DARCOS < C 1 /C2 ) 

HO=Cl/( 4,0D0 *F*USIN( ANGLO ) 

XO-HO/DS IN ( BETA ) 

WRITE ( 67 400 ) HO»XO 

TRANSFER NORMALIZED DISTANCE TO ACTUAL DISTANCE 

DISTX1-XN1 KXO 
D 1 3TX2 - XN2 KXO 

CALCULATE RANGE DX 

DX-DISTX2-DISTX1 



CALCULATE CONSTANTS FOR SUBROUTINE 

F’AI =4 . ODO^DATAN ( 1 ,ODO) 

WN~2 . ODO^F’Al *F/C1 
C2 1 -C2/C1 

RC2.L=C2*R02/(C1*R01 ) 

ANGLO = DF_ 1 A-DARS 1 N ( H/X ) 

WRITE <6>200) 

CALCULATE PRESSURES WITH RESF'ECT TO DISTANCE XX, 
WHERE XX IS CHANGED DY INCREMENT 



L~N F 1 

DO 101-1 fL 

XX ( I ) =D I S TX 1 \ DFLOAT ( I -1 ) *DX/DFLOA T ( N ) 

D 1 S TX : -XX ( I )/X 

CALL PRES ( DISTXfPAMPfPHAS) 

PA< I )=PAMP 
PH ( I ) - F’HAS 

WRITE ( 6 f 300 ) XX < I ) f PA ( I ) fPH ( I ) 

10 CONTINUE 

WRITE (6^800) 

PLOT THE PRESSURE AMPLITUDE US, DISTANCE XX 
AND PRESSURE PHASE US, DISTANCE XX 

WRITE ( 67 700) Cl fC2fR01 fR02fXN1 f XN2f X f H» BETA r F 
CALL DPLTP (XXfPAfNfO) 

WRITE ( 678 OO) 

WRITE (6 7 700) Cl 7C2 fR01 7 RG2 r XN 1 7 XN2 7 X 7 H 7 BETA 7 F 
CALL DPLTP ( XX 7 F‘H 7 N 7 0 ) 
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CJ u 



I/O FORMATS 



C 



ICO FORMAT < 5F 1 5*0) 

200 FORMAT ( //22X ? * DISTANCE ' f 17X f ' PRESSURE AMP . ' f 12X f ' PHASE ANGLE'//) 
300 FORMAT ( 20X f D 1 5 * 7 f i OX > 0 i 5 » 7 f 1 OX f D 1 5 ♦ 7 ) 

400 FORMAT ( //20X ? ' LUUES T POSSIBLE MODE CUT OFF DEPTH ' r ?X f D L 5 . 7 
C /20Xf 'DISTANCE FROM APEX ' ? 20X 7015.7//) 



GOO FORMAT (SF15.0) 

600 FORMAT ( C 3 ) 

/OO FORMAT ( //25X t ' INPUT 
C 25X > ' Cl 

C 25X t ' RO 1 = 

C 25 X f ' XN1 = 

C 25 X f ' X =■• 

C 25X7 'BETA* 

GOO FORMAT (1H1) 

STOP 

END 



DATA' f 1 5X t ' UFDGEO ' // 
f D15.7 f5Xf ' C 2 *' r 015*7/ 



f D 1 5 ♦ 7 f5X > 'R02 
? D 1 5 « 7 f 5X f ' XN2 
f D15 . 7 f 5X f 'H 
7 0 15 • 7 7 5X 7 ' F 



*' fLU 3 . 7 / 

= ' f D 1 5 ♦ 7 / 
*'7015.7/ 
*' 7015.7/) 
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c 

c 

C *** SUB PROGRAM *** 

C 

c * * * * * * * * * * .* * * * *# * * # * * * * * * * * * * * * * :# K it * * * * * .* * * 1c * * * * * * * * * * * * * * * t * * :* * It * * * * * 
C * SUBROUTINE 'PRES' CALCULATES THE PRESSURE AMPLITUDE AND PHASE AT A* 
C *POINT ON THE BOTTOM BY USING IMAGE THEORY . M 

C 

C WHERE C21*RC21*BETAfUN*PAI*X*H AND F ARE INPUT CONSTANTS 
C FROM MAIN ROUTINE* 

C DISTX IS INPUT VARIABLE FROM MAIN ROUTINE t AND r 

C PAMP AND PfIAS ARE OUTPUTS TO MAIN ROUTINE ♦ 

C 

C * NOTATION 
C 

C 



C SYMBOL MEANING , I/O CHARACTER 

C 



c 


PRES 


PRESSURE CAUSED BY N TH PAIR OF IMAGES 


P* Z 


c 


p 


SUM OF PRES . < IF ALL POSSIBLE PAIR OF IMAGES ARE 


P>Z 


c 




ADDED P IS PAMP) 




c 


ANGLN 


ANGLE FORMED BY IMAGE r APEX AND BOTTOM 


P*R 


c 


DD 


PARAMETER FOR CONVENIENCE 


P r R 


c 


R 


TOTAL REFLECTION LOSS OF (N')TH IMAGE 


PfZ 


c 


N 


N TH IMAGE OR N TH PAIR OF IMAGES 


P* I 


c 


M 


M TH BOTTOM BOUNCE 


P* I 


c 


MM 


MAXIMUM NUMBERS OF BOTTOM BOUNCES 


P* I 


c 


THE T AO 


ANGLE FORMED DY THE N TH IMAGE* THE RECEIVING POINT 


P*R 


c 




AND THE BOTTOM 




c 


DR 


DISTANCE FROM N TI-1 IMAGE TO THE RECEIVING POINT 


P*R 


c 




NORMALIZED BY DIVIDING BY XO 




c 


THETAM 


GRAZING ANGLE TO THE BOTTOM 


P * R 


c 




(INTERMEDIATE REFLECTION) 




c 


CHECK 


IDENTIFIES IF GRAZING ANGLE EXCEEDS CRITICAL ANGLE 


P r R 


c 


PSA I 


PARAMETER IN REFLECTION COEFFICIENT 


P*Z 


c 


REFL 


RE- LECTION COEFFICIENT GF INTERMEDIATE BOUNCE 


PfZ 


c 


REFLNO 


REFLECTION COEFFICIENT OF THE LAST BOUNCE 


PfZ 


c 


DIRE T 


SOURCE DIRECTIVITY 


pf R 


c 

c 


Z 


PARAMETER FOR CONVENIENCE 


P*Z 



0 
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non n n n o n noon non non non non noon n m> no o non noon 



* # * * * # :* * * * * 1 # * * * * * % % * * * * * * * * .* * * * * 4c ¥ * * * * ft # * * A * Y M ****************** * **** 



< (PROGRAM) ) 

SUBROUTINE PRES ( DISTX fPAMPfPHAS) 

DEFINE CHARACTERS OF VMKI ABLE f PARAMETERS AND CONSTANTS 

CMPLICIT REAL*U< A-MfO-Z) > INTEGER ( I -N ) 

COMPLEX * 1 6 PSA I , REF- L f REFLNO rfitZ, PRES f P 
COMMON C2 1 t RC2 1 f BETA f WN f PAI f ANGLO f X f H f F 

INITIALIZE COMPLEX PRESSURE P=0.0+J0.0 

P-DCMPLXCO .ODOfO .ODO) 

RESET COUNTER N=0 



N=0 

10 CONTINUE 

CALCULATE ANCLN AND DETERMINE IF ANGLE IS IN THE RANGE OR NOT 
IF ANGLE. GT. PA I GO TO THE NEXT STEP 

ANGLN-2 . ODO Y I NT ( (N+l * )/2* )*BETA+<-l ) YYNYANGLO 
IF (ANGLN.GE.PAI) GO TO SO 

CALCULATE PARAMETER THETAO AND DR 

El 0 s DC OS ( ANGLN ) -DI S TX 
THETAO -DA T AN? ( DS IN ( ANGLN > f DO > 

DR=DSGR TCt. ODO f D 1 S TX**2-2 . 0D0*0 ISTX*DCOS ( ANGLN ) ) 



INITIALIZE THE TOTAL REFLECTION COEFFC IENT R=1.0+J0.0 
R=DCMPLX< 1 .ODOfO* ODO ) 

DEFINE MM BY EACH PATH WHICH IS DETERMINED DY THE N TH IMAGE 

MM=IDIN T< ANGLN/ <2 .ODO*DETA> ) 

L =MM+ 1 
DO 40 M- 1 t L 

CALCULATE PARAMETERS THE TAM AND CHECK 

THETAM-rHETA0-2.0D0*DFL0AT(M-l ) YBETA 
CHECK- 1 ♦ 0 D 0- C 2 1 * * 2 Y D C OS ( T HE T A M ) A * 2 

IDENTIFY THE INCIDENT ANGLE WHICH IS LESS THAN THE CRITICAL 
ANGLE OR NOT 

IF ( CHECK. GT.O. ODO) GO TO 20 

CALCULATE PARAMETER PSAI FOR THE REFLECTION COEFFC IENT OF EACH 
BOUNCE. THERE ARE TWO WAYS WHICH DEPEND ON THE IDENTIFICATION 
OF A CHECK 

PSAI-DCMPLX < 0 . ODO f -DSGR T ( -CHECK ) /DS IN ( THETAM ) ) 

GO TO 30 
20 CONTINUE 

PS AI -DCMPLX < D SORT ( CHECK ) /DS IN ( THETAM ) f 0 • 0 D 0 ) 

30 CONTINUE 

CALCULATE REEL 

REFL = (RC21-PSAI )/( RC21 + PSAI ) 
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non non non n n n n n n n n n non 



IF THE BOUNCE IS THE LAST FROM THE SOURCE t THEN REFLNO=REFL 
I F ( M . E 0 ♦ 1 ) REFLNO=REFL 

CALCULATE THE TOTAL REFLECTION COEFFC I ENT R FOR EACH PATH 



R-R5KREFL 
<10 CONTINUE 

CALCULATE PARAMETERS Z AND PRES 
UMPIRE DIRE r IS THE DIRECTIVITY OF THE SOURCE 
DIRET IS USED ONLY WHEN THE SOURCE IS DIRECTIONAL 
OTHER WISE Ei I RET = 1 * 0 

Z-DCMPLX ( 0 ♦ ODO r -WN KX*DR ) 

0 IRET~DS I N ( 0 » 0 1 ODO*UN'KDSIN ( IHETAM- ANGLO ) )/ 

C ( 0 . 01 ODO *UN'f QSIN( THETAM- ANGLO) ) 

PRES=BIRET*CDEXP(Z)*<-1 ) **INT< (N+l . )/2. >* 

C ( 1 • 0 0 0 t- 1 .ODO/REFLNO MR/DR 

CALCULATE THE PRESSURE P WHICH IS A SUM OF PRES 



P-P+PRES 



SET TO THE NEXT PAIR OF IMAGE 

N = N+1 
GO TO 10 
50 CONTINUE 

CALCULATE PA AND PH f AND RETURN BACK TO THE MAIN ROUTINE 



FAMP=CDAB?j<P> 

PREAL=P 

P-P-PREAL 

PP~P>KDCMPLX ( 0 » ODO 9 1 .ODO) 
PI MACROS IGN< CDABS(P) f-PP) 
PHAS = DA T AN2 ( P L MAG > PREAL ) 
RETURN 
END 
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APPENDIX B 



INCIDENT ANGLE FOR A DISTANCE SOURCE 

General expression of an incident angle is given 
by Eg. (II 39-a) 



0 = 9 2mS 

nm no 



(II 39-a) 



where the maximum number of bounces M is limited by 






M = INT [j|] 



(II-40; 



For a distance source 



0 ~ <p 

mo n 



(11-48) 



A few incident angles 9 for a distance source are shown 

nm 

in Table B-l. From this table it is seen that the reflec- 

til 

tion coefficient for the 4^ image is 



r ( 4 ) = r(0 4 Q )r(0 41 )r( 0 42 ) = r (2) r(0 4Q ) 

and by mathematical induction Eqs. 11-50 can be deduced. 
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m 

n 


0 


1 


2 


3 


• • ♦ 


0 












1 


28-^0 










2 


2B+<J> 0 










3 


4B-4> 0 


2 ' 8-*o 








4 


4 6+4> 0 


23 + 4> 0 








5 


66-<D 0 


4 8-*o 


2B -*o 






6 


6S+4> 0 


4 8+<t> 0 


23+^0 






7 


«"»o 


63-^0 


4B-4, 0 


2B-$ o 




• 













Table B-l 0 

run 
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